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Abstract 

Pulmonary adenoma susceptibility 1 {Pasl) is the major locus responsible for lung tumor susceptibility in mice; among the 
six genes mapping in this locus, Kras is considered the best candidate for Pasl function although how it determines tumor 
susceptibility remains unknown. In an (A/J xC57BL/6)F4 intercross population treated with urethane to induce lung tumors, 
Pasl not only modulated tumor susceptibility (LOD score = 48, 69% of phenotypic variance explained) but also acted, in 
lung tumor tissue, as an expression quantitative trait locus (QTL) for Kras-Ak, one of two alternatively spliced Kras 
transcripts, but not Kras-AB. Additionally, Kras-Ak showed differential allelic expression in lung tumor tissue of (A/J xC57BL/ 
6)F4 heterozygous mice, with significantly higher expression from the A/J-derived allele; these results suggest that c/s-acting 
elements control Kras-Ak expression. In normal lung tissue from untreated mice of the same cross, Kras-Ak levels were also 
highly linked to the Pasl locus (LOD score = 23.2, 62% of phenotypic variance explained) and preferentially generated from 
the A/J-derived allele, indicating that Pasl is an expression QTL in normal lung tissue as well. Overall, the present findings 
shed new light on the genetic mechanism by which Pasl modulates the susceptibility to lung tumorigenesis, through the 
fine control of Kras isoform levels. 
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Introduction 

The Pulmonary adenoma susceptibility 1 (Pasl) locus, mapping 
in the distal region of chromosome 6, is the major modulator of 
lung tumor susceptibility in mice [1,2]. In inbred mouse strains, 
the Pasl locus has a conserved haplotype consisting of six genes 
clustered in a ~450-kb region [3] . From proximal to distal, these 
genes are branched chain aminotransferase 1 (Bcatl), lyrnphoid- 
restricted membrane protein (Lrmp), cancer susceptibility candi- 
date 1 (Cascl), LYR motif containing 5 (Lyrm5), »-Ki-ras2 Kirsten 
rat sarcoma viral oncogene homolog (Kras) and intermediate 
filament tail domain containing 1 (Ifltdl). Four of these genes, 
namely Lrmp, Cascl, Kras and Ifltdl, have been singled out as 
candidate genes for Pasl locus functions in fine-mapping studies 
[4-6]. Although for none of these genes is there a clear 
demonstration of a role in lung tumor susceptibility, there is 
substantial evidence supporting the involvement of Kras as the key 
effector of the Pasl locus. 

Kras encodes a small GTPase that functions as a molecular 
switch in signal transduction, influencing cell proliferation [7]. 
Permanendy activating mutations, at codons 12, 13 and 61, are 
frequently found in both spontaneous and chemically induced 



lung tumors in mice [8,9] . Moreover, in mouse models in which 
mutant Kras can be activated by somatic recombination in the 
lung, animals are highly susceptible to lung tumorigenesis and 
develop multiple lung tumors at 100% incidence with a short 
latency [10,1 1]. Nonetheless, as heterozygous Kras knockout mice 
have higher susceptibility to chemically induced lung tumorigen- 
esis than wild-type mice, the wild-type Kras allele may have a 
tumor suppression function [12]. The double role of Kras in lung 
carcinogenesis — tumor suppressor when wild-type and oncogene 
when mutated — led to the hypothesis that lung cancer suscep- 
tibility could result from the subtle balance between expression 
levels of wild-type and mutated Kras [13]. Given the frequent 
occurrence of activating Kras mutations in mouse lung tumors, an 
increase in Kras gene expression in lung tumor tissue could be 
expected to raise the level of active Kras protein, thereby 
providing the growth advantage characteristic of neoplasms. 
Evidence supporting this mechanism was provided by the 
observation that, in mice, Kras mRNA levels in normal lung 
tissue were ~2-fold higher in strains susceptible to lung 
tumorigenesis (both highly susceptible A/J mice and intermedi- 
ate-susceptible FVB/N and 129/Sv mice) than in a resistant 
strain (C57T3L/6) [13]. 
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Author Summary 

A person's risk of developing cancer depends on both 
genetic and environmental factors. To study the genetic 
predisposition to cancer without the influence of environ- 
mental variables, scientists study mice treated with 
urethane, a chemical carcinogen that induces lung tumors. 
By crossing inbred (genetically identical) strains of mice 
that are either resistant or susceptible to urethane-induced 
cancer, researchers can search for genes associated with 
tumor formation in the offspring. From previous work of 
this type using second-generation mice, it was already 
known that a region on chromosome 6 was associated 
with tumor formation. Now, a new study, carried out in a 
fourth-generation mouse population, focused to a single 
gene of chromosome 6 called Kras. This gene forms two 
different messenger RNA transcripts, called Kras-AA and 
Kras-AB, that produce two proteins with slightly different 
structure and, perhaps, function. The study found that 
mice susceptible to lung tumors have relatively more Kras- 
4A messenger RNA than resistant mice and that this 
difference may be due to small variations in the DNA near 
or within this gene. 

Our understanding of the mode of action of Kras in lung 
tumorigenesis is further complicated by the existence of two main 
transcripts, namely Ai&s-4A and -4B, generated by alternative 
splicing of its fourth coding exon. The two transcripts differ in 
their 3 '-termini and give rise to two proteins (of 189 and 188 
residues, respectively) with different C-terminal sequences. Be- 
cause the C-termini of Ras proteins function in plasma membrane 
binding (through both electrostatic binding of basic residues and 
hydrophobic binding of fatty acylated residues) [14], sequence 
variations in this region could affect the biological functions of the 
isoforms. Indeed, in transfected cells, human Ara.r-4A, but not Kras- 
4B, efficiently induced transformed foci and enabled anchorage- 
independent growth [15,16]. Additional differences between the 
two isoforms regard their expression levels, which were found to be 
modulated during mouse embryogenesis and to vary in different 
tissues [17,18]. Moreover, the ratio between Kras-4A and Aras-4B 
mRNA levels in lung tissue was higher in mouse strains susceptible 
to lung cancer than in resistant ones [17]. More recent evidence 
implicating the Kras-4A oncoprotein in lung tumorigenesis was 
provided by two studies in which mice expressing only the Kras-iB 
isoform had greater resistance to chemical carcinogenesis than did 
wild-type mice also expressing Kras-4A [19,20]. 

Despite these observations, several aspects of the candidacy of 
Kras as the major effector of the Pasl locus in lung tumor 
susceptibility remain to be clarified. Wild-type alleles from 
susceptible and resistant mice code for the same protein; 
consequendy, the mechanism by which Kras determines genetic 
susceptibility to lung tumorigenesis could, instead, depend on its 
overall expression level or on the ratio of its isoforms in normal 
lung tissue. 

To address these points and to clarify the involvement of the 
two Kras isoforms, we studied Kras mRNA levels and germline 
variants using an established model of urethane-induced lung 
tumorigenesis in an advanced intercross population between A/J 
mice (tumor susceptible) and C57BL/6 mice (resistant). This 
model was previously instrumental in defining the Pasl locus based 
on the pattern of tumorigenesis in F2 intercross mice [1]. In the 
present study, we used F4 mice (of a new pedigree) in order to 
have a greater resolution power. After confirming Pasl as a 
quantitative trait locus (QTL) modulating lung tumor multiplicity 



in this pedigree, we performed a genome-wide association study to 
look for expression QTL modulating the levels of Kras-iA and 
Kras-4B transcripts in lung tumors and normal lung tissue. Our 
findings indicate that KrasAA and, to some extent, also Ara.s-4B are 
modulated by m-acting elements within the Pasl locus itself. These 
results underpin a possible causal relationship between germline 
variations and mRNA levels of the Ara.r-4A isoform, consequendy 
influencing lung tumor susceptibility. 

Results 

A/J and C57BL/6 mice were mated and bred to the fourth 
generation (ABF4) to establish a pedigree for the current study. 
Male ABF4 mice (n = 183) were treated with a single injection of 
urethane at 4 weeks of age to induce lung tumor formation, 
sacrificed 36 weeks later and assessed for lung tumor multiplicity 
(Nlung). Values of Nlung ranged from 0 (in 27 animals) to a 
maximum of 33. 

To assess the genetic control of susceptibility to lung tumori- 
genesis in this pedigree, we performed genome-wide linkage 
analysis on the 183 urethane-treated ABF4 mice using Illumina 
SNP-arrays, which permit the genotyping of up to 1449 single 
nucleotide polymorphisms (SNPs). After quality control filtering, 
genotype data were obtained for 548 informative (polymorphic) 
non-redundant SNPs dispersed over the whole genome. Simple 
interval mapping was performed to detect QTLs associated with 
squared root-transformed values of Nlung. This analysis identified 
a single major locus in the telomeric region of chromosome 6 
(Figure 1A). The QTL peak had a LOD score of 48.2 and was 
centered around marker rs6265387 [Pthlh gene region) near the 
Pasl locus; moreover, at the Kras 37-bp marker, located 2-Mb 
proximal to rs6265387 and distinguished from the SNP by 3 
recombination events, the LOD score was almost identical 
(LOD = 48.0). The QTL accounted for 69% of the total 
phenotypic variance in the ABF4 population. When the mice 
were grouped according to genotype (Figure IB), Nlung in the 59 
animals homozygous for the A/J-derived allele (represented by a 
G at rs6265387) was 1.5-fold higher than that in the 82 
heterozygous mice and ~ 24-fold higher than that in the 42 mice 
homozygous for the C57BL/6-derived allele (represented by an A 
at rs6265387; P< 1.0 x 10" 6 , ANOVA followed by Tukey's test for 
multiple comparisons). 

These results, which confirm the key role of the Pasl locus in 
murine lung tumorigenesis, establish this new intercross 
population as suitable for studying Pa.fi function. In this cross, 
lung carcinogenesis can be considered a monogenic trait due to 
the overwhelming genetic effect of the Pasl locus. Therefore, 
we used this pedigree to test whether the Pasl locus exerts its 
effect on lung tumorigenesis through a modulation of Kras 
mRNA levels. To this aim, we performed linkage analyses to 
identify expression QTLs associating with the levels of Kras-4A 
and KrasAB mRNA. 

Pasl is an expression QTL controlling Kras-AA transcript 
levels in lung tumors of urethane-treated mice 

From 80 of the 183 ABF4 urethane-treated mice used in the 
genetic linkage analysis, we were able to resect a tumor specimen 
from lung tissue. This subgroup included 37 mice homozygous at 
rs6265387 for the A/J-derived allele (associated with tumor 
susceptibility), 34 heterozygous mice, and 9 mice homozygous for 
the C57BL/6 resistant allele that nevertheless had developed a 
lung tumor. Indeed, the intrinsic low susceptibility to lung 
tumorigenesis of the 42 animals homozygous for the C57BL/6- 
derived allele (only 9 of these mice developed a lung tumor) did 
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Figure 1. Lung tumor multiplicity in 183 urethane-treated male ABF4 mice is controlled by the Pasl locus, (a) LOD score plot for 
chromosome 6 on which a quantitative trait locus (QTL) for lung tumor multiplicity (square root transformed values) mapped to the telomeric region. 
The QTL peak (LOD score = 48, phenotypic variance explained = 69%) overlapped with the Pas! locus. Tick marks show the position of 37 genotyped 
markers, including rs6265387 at the QTL peak. Horizontal line indicates the 95% LOD threshold, (b) Number of lung tumors per animal, grouped 
according to genotype at rs6265387. Mice homozygous for the A/J-derived allele (GG; n = 59) had more tumors than either heterozygous animals 
(n = 82) or mice homozygous for the C57BL/6-derived allele (AA; n = 42). ***P<1.0x1(T 6 versus the A/J-derived allele, ANOVA followed by Tukey's 
test for multiple comparisons. The line within each box represents the median; upper and lower edges of each box are 75 th and 25 th percentiles, 
respectively; upper and lower bars indicate the highest and lowest values less than one interquartile range from the extremes of the box. 
doi:10.1371/journal.pgen.1004307.g001 



not allow us to analyze a subgroup of the original population with 
a genotype ratio typical of such intercrosses (i.e., 1:2:1). 

RNA was extracted from each specimen and used in 
quantitative PCR to measure the levels of KrasAA and -4B 
mRNA. Genome-wide linkage analysis for KrasAA found a single 
expression QTL on chromosome 6 (peak LOD = 4.5), corre- 
sponding to the Pasl locus (Figure 2A). In contrast, no expression 
QTL was found for KrasAB on any chromosome, including 
chromosome 6 where LOD scores remained below the threshold 
for significance. When the 80 mice were grouped according to the 
genotype at rs6265387, we found that KrasAA mRNA was highest 
in mice homozygous for the A/J-derived susceptible allele (GG), 
intermediate in heterozygous (AG) animals, and lowest in mice 
homozygous for the C57BL/6 resistant allele (AA); these 
differences were significant (GG vs. AG, P= 3.2x10 ' 5 ; GG vs. 
AA, P= 8.9 x 10" 3 ; ANOVA followed by Tukey's test for multiple 
comparisons, Figure 2B). These results implicate KrasAA in the 
mechanism of Payi-dependent tumor formation. 

Pas 7-dependent tumor susceptibility correlates with 
Kras-4A transcript levels in normal lung tissue 

The genetic susceptibility to cancer is an intrinsic feature of 
normal tissue that, in urethane-treated mice, influences both the 
probability of tumor initiation and the number of tumors that 
develop. Consequendy, it is likely that the fundamental elements 
underlying this phenomenon are present not only in tumors but 
also in normal tissue. To test this hypothesis, we attempted to 
validate in normal lung tissue the genetic linkage we observed in 
lung tumors. Therefore, we genotyped 1 1 1 untreated male ABF4 
mice for nine markers on chromosome 6 spanning from 96.7 Mb 
to 148.3 Mb. These markers include 8 SNPs and one 37-bp 



sequence variation that in C57BL/6 mice occurs as a tandem 
repeat; this insertion mutation is located in the second intron of 
Kras [21]. Additionally, we assayed lung mRNA from these mice 
for the two Kras isoforms, and examined the association of mRNA 
levels with genotype. 

A strong linkage was found between the level of KrasAA mRNA 
and genotyped markers in the Pasl locus, describing a LOD curve 
with a peak of 23.2 (Figure 3A); this expression QTL explained 
62% of the phenotypic variance. A much weaker, yet statistically 
significant linkage for the KrasAB transcript was observed 
(maximum LOD score = 4. 1). These results confirm and strength- 
en the major role of KrasAA isoform found in tumors (Figure 2). 
We then grouped the mice according to genotype at the Kras 37-bp 
insertion to examine the relationship between Kras mRNA levels 
and the number of inherited A/J-derived alleles of the Pasl locus. 
A strong association was found for KrasAA levels (P= 9.7 x 10 14 , 
ANOVA; square-root-transformed values are shown in 
Figure 3B): mice with two A/J-derived alleles (i.e. no insertion, 
-/-) had on average 1.9-times more transcript than did mice 
receiving two C57BL/6-derived alleles (ins/ins), and heterozygous 
mice had intermediate levels. In contrast, the linkage between Pasl 
genotype and KrasAB transcript levels was weaker (P= 9.7 x 10~ 5 , 
ANOVA, Figure 3C). 

We then examined whether this specific genetic control exerted 
by the Pasl locus was also present in normal lung of the parental 
inbred strains A/J and C57BL/6 (Figure 4). The median level of 
the KrasAA transcript in A/J (tumor-susceptible) mice was 2.3-fold 
higher than that of C57BL/6 (resistant) mice (P= 6.0 x 10~ 9 , 
ANOVA; square-root-transformed values are shown in 
Figure 4A). In contrast, the KrasAB levels were not significantly 
different between mouse strains (Figure 4B). 
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Map position (cM) rs6265387 

Figure 2. Kras-4I\ levels in lung tumors from 80 urethane-treated ABF4 mice are controlled by Pas 1 locus, (a) Expression QTL analysis of 
Kras transcripts showed that square-root-transformed levels of Kras-Ak linked to the Pasl locus (LOD = 4.5). No significant linkage was observed for 
the Kras-4B mRNA isoform. Tick marks show the position of the genotyped markers in a recombinational map. Horizontal line at LOD = 2.55 marks the 
95% threshold for significance, (b) Relative expression levels of the Kras-Ak isoform according to rs6265387 genotype. Mice homozygous for the A/J- 
derived susceptible allele (GG, n = 37) had higher levels than either heterozygous animals (AG, n = 34) or mice homozygous for the C57BL/6 resistant 
allele (AA, n = 9). ***P<0.001, **P<0.01 vs. GG mice, ANOVA followed by Tukey's test for multiple comparisons. Values are means and SE. 
doi:1 0.1 371 /journal.pgen.1 004307.g002 



Altogether, these results indicate that the A/J-derived allele of 
the Pasl locus, which confers susceptibility to lung tumorigenesis, 
is also associated with higher steady-state levels of the KrasAA 
splice variant. These observations suggest that subtle modulation 
of KrasAA mRNA production or stability may be the key effector 
of Pflii-controlled lung tumor susceptibility, possibly via a ex- 
acting element within Pasl itself. 

Differential allelic expression of Kras-Ah is due to 
sequence variants in c/s-regulatory elements 

To test whether the allele-specific levels of Kras isoforms are 
attributable to variations in ciy-regulatory elements also mapping 
in the Pasl locus, we examined the possibility of differential allelic 
expression. For this analysis, we considered two SNPs (rs29968550 
and rs30022167, mapping in a region of Kras mRNA common to 
the two isoforms) whose genotypes permitted us to distinguish the 
A/J- from the C57BL/6J-derived allele. By pyrosequencing, we 
determined the frequencies of the two alleles in Kras-4A and -4B 
cDNA from normal lung tissue of 20 heterozygous untreated 
ABF4 mice and from lung tumor specimens of 15 heterozygous 
treated mice and in genomic DNA from the same animals, used as 
control for the amplification efficiency of the two alleles (Figure 
SI). Indeed, when allelic imbalance is tested, uneven amplification 
unrelated to differential transcription may occur in some of the 
assays. Therefore, the experimentally obtained allelic ratio for 
genomic DNA, which in principle is equal to 1 in heterozygotes, is 
used as the baseline to evaluate the ratios observed for cDNAs. 

Based on the results for both SNPs, comparing the mean allelic 
ratios for Kras-4:A cDNA with those of the corresponding genomic 
DNA, we observed a ~2-fold higher expression of the A/J-derived 
allele in Kras-4A cDNA in both tumor and normal tissue indicating 



a higher steady state of this isoform associated with the A/J- 
derived allele, possibly due to either a higher production rate or a 
lower degradation rate (Figure 5). For Ara.r-4B, the mean allelic 
ratios for cDNA from tumors were similar to those for genomic 
DNA, meaning equal production from both alleles; in normal 
tissue, however, a slightly higher ratio for Ara,r-4B cDNA relative to 
genomic DNA was observed (Figure 5). 

In the statistical analysis, allelic ratios were log m-transformed to 
approximate a normal distribution before applying Welch's t test 
(Table 1). This analysis showed that the high allelic ratios for Kras- 
4A cDNA, observed at two SNPs and in both normal and tumoral 
lung tissue, are significantly different from those of genomic DNA 
(P<0.001, see Table 1). In the case of KrasAB cDNA, the allelic 
ratios were significantly different from that of genomic DNA only 
for normal lung tissue (P<0.001). These results provide evidence 
for the differential allelic expression of Kras transcripts, with higher 
levels being produced from the A/J-derived allele in ABF4 
intercross mice. These findings are compatible with the presence 
of allele-specific germline variations in m-acting elements that 
mainly influence the selection or stability of the KrasAA isoform. 

Discussion 

In this study, we verified the major role of the Pasl locus in 
modulating lung tumorigenesis in an ABF4 advanced intercross 
population, where the lung tumor multiplicity (Nlung) phenotype 
behaved as a monogenic trait under the control of the Pasl locus, 
explaining almost 70% of the phenotypic variance. In addition, we 
found that Pasl is an expression QTL because it controls the level 
of the Kras-4:A isoform in urethane-induced lung tumors. This 
genetic control is an inherited trait, in that it was already present in 
normal lung tissue of untreated mice. Finally, we observed the 
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A B C 




Genotyped SNPs Kras-37bp Kras-37bp 



Figure 3. Kras-4/k expression is highly controlled by Pasl locus in normal lung, (a) Genetic linkage analysis of expression levels (square root 
transformed) of two Kras transcripts in 1 1 1 untreated ABF4 mice showed that the Kras-4A isoform was strongly linked to the Past locus (LOD 
score = 23.2) whereas the KrasAB isoform showed a weaker linkage (LOD score = 4.1). Horizontal line at LOD = 2.1 3 marks the threshold for 
significance. Tick marks show the position, in a recombinational map, of the nine genotyped markers spanning from chromosome 6 position 96.7 Mb 
to 148.3 Mb. (b, c) Expression levels of Kras-4A and -4B in normal lung tissue, by genotype for the 37-bp insertion mutation common to both 
isoforms. The A/J-derived susceptible allele is negative (-) for the insertion whereas the C57BL/6 allele is positive (ins); 30 mice were -/-, 52 ins/-, and 
29 ins/ins. (b) For Kras-Ak, P = 9.7x1CT 14 , ANOVA. Tukey's test for multiple comparisons, *P = 0.01, ***P = 7.4x1 (T 14 vs. -/-. (c) For Kras-4B, 
P = 9.7x10~ 5 , ANOVA. Tukey's test for multiple comparisons, **P= 1.5x1CT 4 vs. -/-. Values are means and SE. 
doi:1 0.1 371 /journal.pgen.1 004307.g003 





A/J 



C57BL/6 



A/J 



C57BL/6 



Figure 4. Kras-4A is expressed at higher levels in susceptible (A/J) than resistant (C57BL/6) strains. Expression levels of Kras-4A (a) and 
Kras-4B (b) isoforms were measured by qPCR in normal lung of 11 A/J and 11 C57BL/6 mice. * P = 7.1 x10~ 5 , ANOVA. The line within each box 
represents the median of square root transformed values; upper and lower edges of each box are 75 th and 25 th percentiles, respectively; upper and 
lower bars indicate the highest and lowest values less than one interquartile range from the extremes of the box. 
doi:1 0.1 371 /journal.pgen.1 004307.g004 
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Tumor Normal Tumor Normal 

rs30022167 rs29968550 

Figure 5. Differential allelic expression of Kras-4I\ isoform indicates the existence of functional polymorphisms in Kras gene. Allelic 
ratios were determined by pyrosequencing for rs30022167 and rs29968550, which map in a region of the Kras 3'-UTR common to both isoforms, on 
genomic DNA and cDNA from normal lung tissue (n = 20) and lung tumor specimens (n = 1 5) from ABF4 heterozygous mice. Values are mean and SE; 
*** P<0.001, two-sided Welch's t test (see Table 1 for complete data). 
doi:1 0.1 371 /journal.pgen.1 004307.g005 



Table 1. Differential allelic expression of Kras isoforms in 
mice. 


lung tumor specimens and normal lung tissues from ABF4 intercross 




Mean log 10 -transformed allelic 
ratio* 


SNP Tissue 


Isoform 


cDNA 


Genomic DNA 


Difference (95% CI) 


P value* 


rs30022167 Tumor 


Kras-4A 


0.26 


-0.05 


0.31 (0.20 to 0.43) 


3.7x10~ 5 




Kras-4B 


0.01 


-0.05 


0.06 (-0.03 to 0.16) 


0.17 


Normal 


Kras-4A 


0.22 


-0.07 


0.29 (0.26 to 0.32) 


<2.2x10~ 16 




Kras-4B 


0.02 


-0.07 


0.09 (0.06 to 0.12) 


6.1 x10~ s 


rs29968550 Tumor 


Kras-4A 


0.42 


0.16 


0.26 (0.13 to 0.39) 


8.4 x10~ 4 




Kras-4B 


0.22 


0.16 


0.06 (-0.03 to 0.15) 


0.19 


Normal 


Kras-4A 


0.37 


0.04 


0.33 (0.30 to 0.37) 


<2.2x10~ 16 




Kras-4B 


0.21 


0.04 


0.17 (0.15 to 0.20) 


<2.2x10~ 16 





* Mean of log 10 transformed allelic ratios (A/J-derived allele/C57BL/6-derived allele) for 20 samples (normal tissue) or 15 samples (tumor tissue). 
^Calculated on differences between cDNA and DNA allelic ratios using a two-tailed Welch's f test. 
doi:1 0.1 371 /joumal.pgen.1 004307.t001 



PLOS Genetics | www.plosgenetics.org 



6 



April 2014 | Volume 10 | Issue 4 | e1 004307 



Pas/ Locus Modulates Kras-Ak Levels 



preferential expression of Kras-4A from the A/J-derived allele in 
heterozygous mice, suggesting the existence of allele-specific 
germline variations in as-acting elements that influence splicing 
bias, selection or stability of this isoform. 

We carried out this genetic study in an ABF4 intercross since it 
was an effective way to accumulate, in a population of a relatively 
small size, 3-fold more recombination events than those possible 
in a conventional F2 intercross population of the same size [22]. 
Indeed, each ABF4 mouse underwent six informative meioses 
compared to the two informative meioses of an F2 mouse, thus 
providing improved resolution for mapping loci affecting strain- 
related phenotypes [22]. This methodological choice allowed us 
to obtain a sharp peak of linkage between Pasl and Nlung, with 
very high association (LOD score = 48 in 183 ABF4 mice). The 
high rate of recombination in this intercross was also useful for 
testing our hypothesis that Kras isoform expression is modulated 
by Pasl. 

Animals homozygous for the C57BL/6-derived allele had low 
susceptibility to lung tumorigenesis and, therefore, rarely 
developed lung tumors after urethane treatment. Hence, the 
genetic class with the resistant phenotype was under-represented 
in the experiment designed to identify expression QTLs in 
tumors. Another difficulty that presented in our research was 
the known dysregulation of gene expression architecture in 
cancer, including mouse lung tumors [23,24], which could have 
reduced the detectability of expression QTLs. Despite these 
limitations, we observed a significant linkage of the Pasl locus 
with mRNA levels of Kras-4A (LOD score = 4.5 in 80 mice) but 
not of Kras-4B. 

In normal lung tissues from untreated mice, where all three 
genotypes were present at roughly the expected ratio and where 
no pathologic alterations could have played a confounding role, 
we observed that the Kras-4A isoform levels were again signifi- 
cantly linked with the Pasl locus (LOD score = 23.2 in 1 1 1 mice). 
In these untreated animals, the Kras-4B isoform also had 
significant linkage with Pasl, although the LOD score of 4. 1 was 
only ~17% that obtained for Kras-4A. 

The differential allelic expression of Kras-4A and, to some 
extent, Kras-4H, in normal lung tissue and in lung tumor 
specimens, indicated the existence of regulatory polymorphisms 
located within the Pasl locus. These cw-acting elements have not 
yet been identified. However, since the Kras gene promoter has 
been shown to have similar activity in susceptible and resistant 
mice [25], it is conceivable that these m-acting elements are 
located in other regions such as introns or the 3 '-untranslated 
region of the Kras transcript. Regulatory elements in these 
regions can be expected to impact more on mRNA splicing and 
stability than on promoter activity. Several polymorphisms in 
non-coding regions of Kras have already been found in various 
mouse inbred strains and the functional roles of some of them 
have been investigated [21,25,26]; these studies, however, 
reached contrasting conclusions. Future investigations could 
focus on the targeted re-sequencing of the A/J-derived Kras gene 
and flanking regions, to search for not yet identified polymor- 
phisms that affect splicing selection or mRNA stability of Kras 
isoforms. The candidate functional polymorphisms could then 
be tested by cloning into reporter vectors for transfection and 
assay of cell lines. 

Overall, these results shed new light on the role of expression 
QTLs in mouse lung tumorigenesis, since they indicate that the 
genetic control exerted by Pasl mostly affects the expression of the 
Kras-4A transcript. These results point to the 4A isoform as the 
functional element of Pasl, the major locus modulating lung 
tumorigenesis in mice. 



Materials and Methods 

Ethics statement 

All animals received humane care according to the criteria 
oudined in a protocol approved during the meeting board of 
December 21, 2006, by the institutional ethical committee for 
animal use (CESA) at the Fondazione IRCCS Istituto Nazionale 
dei Tumori. 

Mouse crosses, tissues, DNA and RNA 

A pedigree of intercross mice was generated by mating lung 
tumor-susceptible A/J mice (A; 5 females) with lung tumor- 
resistant C57BL/6 mice (B; 5 males). From generation ABF2 to 
ABF4, male and female mice were individually labeled, randomly 
selected using a random number generator and bred, avoiding the 
pairing of siblings (20 female and 20 male ABF2 mice, and 55 
female and 55 male ABF3 mice were mated). After weaning, ABF4 
mice were sexed and ear-tagged, and a section of tail was collected 
for DNA extraction. Male ABF4 mice were then randomly 
assigned to two groups: untreated mice (n = 131) were raised under 
standard conditions, while 188 male mice were treated with a 
single intraperitoneal injection of urethane (1 g/kg body weight) at 
4 weeks of age to induce the development of lung tumors [27]. 

To obtain normal lung tissue, 1 1 1 untreated AFB4 mice, 1 1 A/J 
and 1 1 C57BL/6 mice (all males) were anesthetized and killed at 
16 weeks of age; lung lobes were isolated and frozen in liquid 
nitrogen. To obtain lung tumor tissue, the 183 urethane-treated 
mice were killed at 40 weeks of age; the chest was opened, the 
trachea was dissected to permit needle access, and the lungs were 
filled with 0.5 ml RNALater solution. Then, the lung lobes were 
removed, placed in RNALater, and kept overnight at 4°C degrees. 
The next day, lungs were examined for surface tumors. For each 
mouse, we recorded the total number of tumors (Nlung) and, 
whenever possible, resected one tumor of about 1-1.5 mm under a 
stereo-microscope. These tumor specimens were frozen and stored 
at — 80°C degrees. 

Genomic DNA was extracted from tail samples using the 
DNeasy Blood & Tissue Kit (Qiagen, Valencia, CA, USA) and 
quantified using Picogreen dsDNA Quantitation Kit (Invitrogen, 
Life Technologies, Paisley, UK). Total RNA was extracted from 
normal lung tissue and from lung tumor specimens using RNeasy 
Midi Kit (Qiagen), purified with RNeasy MinElute Cleanup 
(Qiagen), and quantified by spectrophotometry (ND-1000 spec- 
trophotometer, NanoDrop, Wilmington, DE, USA). RNA integ- 
rity was verified using the RNA 6000 Nano Assay Kit (Agilent 
Technologies, Palo Alto, CA, USA); the mean RIN value was 8.9 
(SD = 0.8; range, 6.7 to 9.9). 

Genome-wide SNP and Pasl locus genotyping 

Genomic DNA from urethane-treated ABF4 mice was used for 
genome-wide SNP genotyping carried out with the GoldenGate 
Genotyping Assay according to the manufacturer's protocol 
(Illumina, San Diego, CA, USA), using the Mouse MD Linkage 
Panel representing 1449 mouse loci. 

Additionally, the Pasl locus in untreated ABF4 mice was 
genotyped for nine genetic markers, including eight SNPs 
(rs6349084, rs33863668, rs31000839, rs31005929, rsl3479063, 
rs33893742, rsl3479082, rs3711088) and a 37-bp tandem repeat 
in Kras [21]. Briefly, a genomic fragment surrounding each marker 
was PCR-amplified in reactions containing 30 ng genomic DNA, 
1 U AmpliTaq Gold (Applied Biosystems, Life Technologies), 1 x 
AmpliTaq Gold buffer (Applied Biosystems, Life Technologies), 
1.5 mil MgCl 2 , 200 uM dNTPs and 5 pmol of each of a pair of 
specific primers (Table SI) in a total volume of 25 u.1. To 
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genotype the eight SNPs, PCR products were pyrosequenced on a 
PyroMark Q96 ID system running PyroMark Q96 ID Software 
(Qiagen). To genotype the Kras 37-bp repeat, PCR amplicons were 
analyzed by 3% agarose gel electrophoresis for fragment size. 

Quantitative real-time PCR 

RNA from normal lung (1 ug) or from lung tumor specimens 
(0.5 ug) was used to synthesize cDNA by reverse transcription 
using the Transcriptor First Strand cDNA Synthesis Kit (Roche, 
Basel, Switzerland). Kras-4A and -4B transcript levels in the cDNA 
were measured in quantitative PCR (qPCR) assays using Fast 
SYBR Green PCR Master Mix (Applied Biosystems) and 300 nM 
intron-spanning primers (Table SI). Hypoxanthine phosphor- 
ibosyltransferase 1 (Hprtl, Table SI) was used as reference gene. 
Relative expression levels were calculated using the comparative 
Ct method using one of the synthesized cDNA samples as 
calibrator. 

Detection of differential allelic expression of Kras-AK and 
-4B transcripts 

The allelic expression of Kras-4A and Ara.s-4B isoforms was 
analyzed taking into consideration two SNPs in the 3'-UTR 
region of the Kras gene (Figure SI): rs30022167 (A or C, for the 
A/J- or C57BL/6-derived allele, respectively) and rs29968550 (A 
or G, for the A/J- or C57BL/6-derived allele, respectively). 
Among the animals heterozygous for markers at the Pasl locus, we 
selected 20 mice from the untreated group and 15 from the 
urethane-treated group, from which we already had genomic 
DNA. In addition, from these animals, we also had cDNA from 
normal lung (n = 20) or from tumor specimens (n= 15), having 
been reverse-transcribed for the qPCR experiments. For each Kras 
isoform, we carried out a first amplification step on the cDNA 
samples using a forward primer located in exon 4A (Kras-AA) or in 
the junction between exon 3 and exon 4B (KrasAB) and a common 
reverse primer located in the 3'-UTR region of Kras gene 
downstream of rs30022167 and rs29968550 (Table SI and 
Figure SI). The resulting PCR product as well as a sample of 
genomic DNA from the same animal were PCR-amplified using 
SNP-specific primers. These PCR products were pyrosequenced 
on a PyroMark Q96 ID system running PyroMark Q96 ID 
Software (Qiagen). For each SNP, the proportions of the two 
alleles present in each sample (genomic DNA, KrasAA cDNA and 
KrasAB cDNA from untreated and tumor-bearing mice) were 
determined from pyrogram peak heights, and allelic ratios (A/J- 
derived allele/C57BL/6-derived allele) were calculated. 
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size [28]. 

Differences between genotype groups were analyzed by one- 
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Supporting Information 

Figure SI Experimental design for measuring the allelic 
expression of Kras isoforms. (a) First amplification step on cDNA 
was carried out using a forward primer located in exon 4A (green 
arrow) for the amplification of KrasAA, or located in the junction 
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common reverse primer (red arrow) located in the 3'-UTR region 
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obtained in the previous step and from genomic DNA of the same 
animals, using SNP-specific forward and reverse primers (blue 
arrows for rs29968550 and orange arrows for rs30022167) and a 
sequencing primer (dashed blue arrow for rs29968550 and dashed 
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